Quantum localized modes in capacitively coupled Josephson junctions 
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We consider the quantum dynamics of excitations in a system of two capacitively coupled Joseph- 
son junctions. Quantum breather states are found in the middle of the energy spectrum of the 
confined nonescaping states of the system. They are characterized by a strong excitation of one 
junction. These states perform slow tunneling motion from one junction to the other, while keeping 
their coherent nature. The tunneling time sensitively depends on the initial excitation energy. By 
using an external bias as a control parameter, the tunneling time can be varied with respect to 
the escape time and the experimentally limited coherence time. Thus one can control the flow of 
quantum excitations between the two junctions. 

PACS numbers: 63.20.Pw, 74.50.+r, 63.20.Ry 
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Josephson junctions are the subject of extensive stud- 
ies in quantum information experiments because they 
possess two attractive properties: they are nonlinear 
devices, and also show macroscopic quantum behavior 
[3, 0- El- The dynamics of a biased Josephson junction 
(JJ) is analogous to the dynamics of a particle with a 
mass proportional to the junction capacitance Cj, mov- 
ing on a tilted washboard potential 



U(<p) = -I c — cos ip- I b ip—, 
Zir Zir 



(1) 



which is sketched in FigJTJb. Here <p is the phase differ- 
ence between the macroscopic wave functions in both su- 
perconducting electrodes of the junction, I c is the critical 
current of the junction, and <I>o = h/2e the flux quantum. 
When the energy of the particle is large enough to over- 
come the barrier AU (that depends on the bias current 
lb) it escapes and moves down the potential, switching 
the junction into a resistive state with a nonzero voltage 
proportional to ip. Quantization of the system leads to 
discrete energy levels inside the wells in the potential, 
which are nonequidistant because of the anharmonicity. 
Note that even if there is not enough energy to classically 
overcome the barrier, the particle may perform a quan- 
tum escape and tunnel outside the well, thus switching 
the junction into the resistive state Thus each state 
inside the well is characterized by a bias and a state- 
dependent inverse lifetime, or escape rate. 

Progress on manipulation of quantum JJs includes 
spectroscopic analysis, better isolation schemes, and si- 
multaneous measurement techniques @, 0, 0, H, IS 0, 
and paves the way for using them as JJ qubits in ar- 
rays for experiments on processing quantum information. 
Typically the first two or three quantum levels of one 
junction are used as quantum bits. Since the levels are 
nonequidistant, they can be separately excited by apply- 
ing microwave pulses. 

So far, the studies on JJ qubits focused on low en- 
ergy excitations involving the first few energy levels of 
the junctions. Larger energies in the quantum dynam- 
ics of JJs give rise to new phenomena that can be ob- 
served by using the already developed techniques for 



quantum information experiments. For instance, it was 
suggested that Josephson junctions operating at higher 
en ergies may be used for experiments on quantum chaos 
[H HH3] • Another phenomenon is the excitation of quan- 
tum breathers (QB) [HI, [13, d , which are nearly degen- 
erate many-quanta bound states in anharmonic lattices. 
When such states are excited, the outcome is a spatially 
localized excitation with a very long time to tunnel from 
one lattice site to another. So far the direct observa- 
tion of this kind of excitations evolving in time has not 
been reported. Up to date evidences of QB excitations 
were obtained s_pectroscopically in molecules and solids 



In this work we consider large energy excitations in 
a system of two capacitively coupled JJs [U HH, Ull- 
We study the time evolution of states when initially only 
one of the junctions was excited. In the low energy sec- 
tor such a state will lead to a beating with a beating 
time depending solely on the strength of the capacitive 
coupling. For larger excitation energies the states per- 
form slow tunneling motion, where the tunneling time 
sensitively depends on the initial energy, in contrast to 
the low energy beating time. We calculate the eigen- 
states and the spectrum of the system, and identify quan- 
tum breather states as weakly splitted tunneling pairs of 
states [13, HH ■ These eigenstates appear in the middle of 
the energy spectrum of the system and are characterized 
by correlations between the two junctions - if one of them 
is strongly excited, the other one is not, and vice versa. 
By exciting one of the junctions to a large energy, we 
strongly overlap with QB tunneling states. Consequently 
we trap the excitation on the initially excited junction on 
a time scale which sensitively depends on the amount of 
energy excited, and on the applied bias. We describe how 
this trapping could be experimentally observed in time 
using the nowadays used techniques for manipulating JJ 
qubits. 

The system is sketched in FigfTJa: two JJs are coupled 
by a capacitance C c , and they are biased by the same 
current Ii,. The strength of the coupling due to the ca- 
pacitor is ( = C C /(C C + Cj). The Hamiltonian of the 
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FIG. 1: (a) Circuit diagram for two ideal capacitively cou- 
pled JJs. (b) Sketch of the washboard potential for a single 
current-biased JJ. (c) Sketch of a breather solution in the 
classical dynamics of the system. 



system is 



H=-^ + -^ + U( ( p 1 ) + U( ( p 2 ) + ^P 1 P 2 , (2) 
2m 2m m 



where 



»' = on (^) • 



(3) 



Pi,a = (C c + Cj) 1^1 (^, 2 - C0 2 ,i). (4) 

When the junctions are in the superconducting state, 
they behave like two coupled anharmonic oscillators. The 
plasma frequency is uj p = y / 2irI c /QoCj(l + C)[l~ l 2 ] 1 ^, 
and 7 = Ib/I c is the normalized bias current. The clas- 
sical equations of motion of the system admit breather 
solutions [29l | . which are time periodic, and for which the 
energy is localized predominantly on one of the junctions 
(Figfljc). These orbits are numerically computed with 
high accuracy by using Newton algorithms [3(| HU . At 
the studied energies the classical phase space is mixed, 
and breathers are located inside regular islands, which 
are embedded in a chaotic layer. 

In the quantum case we compute the energy eigen- 
values and the eigenstates of the system. We neglect 
quantum escape for states which will not escape in the 
classical limit. Thus we use a simple model for the single 
J J, where the potential energy is changed by adding a 
hard wall which prevents escape: 



U(lfi) if (f < 7T — ip 
OO if if > 7T — iff) 



(5) 



where ipo — arcsin7 is the position of the minimum of 
the potential and tt — tpo gives the position of the first 
maximum to the right from the equilibrium position ipo 



(Figffl-b). We will later compare the obtained tunneling 
times with the true state dependent escape times. 

The Hamiltonian of the two-junctions system is given 
by H = H! + H 2 + CV, where H t = Pf /2m + U q (<p l ) is the 
single-junction Hamiltonian and V = P\P 2 /m is the in- 
teraction that couples the junctions. The eigenvalues e ni 
and eigenstates \rii) of the single-junction Hamiltonian 
Hi were computed by using the Fourier grid Hamiltonian 
method [13]. Note that |n,) is also an eigenstate of the 
number operator hi with eigenvalue m. In the harmonic 
approximation hi = a\di, where of and hi are the bosonic 
creation and annihilation operators. Since only states 
with energies below the classical escape energy (barrier) 
are taken into account, the computed spectra have a fi- 
nite upper bound. The perturbation V does not conserve 
the total number of quanta n\ + n 2 , as seen from the de- 
pendence of the momentum operators on the bosonic cre- 
ation and annihilation operators in the harmonic approx- 
imation: P1.2 = (<3? /27t)^(1 + QCjfuv p /2(a 1 ^-a\ 2 )/i. 

The Hamiltonian matrix is written in the basis of prod- 
uct states of the single-junction problem {\nx,n 2 ) = 
® |^2)}- The invariance of the Hamiltonian under 
permutation of the junction labels allows us to use sym- 
metric and antisymmetric basis states 

\n>i,n 2 )s,A = —^(\ni,n 2 ) ± |n 2 ,ni)) (6) 

to reduce the full Hamiltonian matrix to two smaller sym- 
metric and antisymmetric decompositions of H , which af- 
ter diagonalization respectively give the symmetric and 
antisymmetric eigenstates of the system. 

In order to identify quantum breather states, whose 
corresponding classical orbits are characterized by energy 
localization, we define the correlation functions: 



/ M (l,2) = (n 1 n 2 ), 
/„(!,!) = (n?)„, 



(7) 

(8) 



where (A)^ = (xp\A\xn), {1x^)1 bein S the set of eigen- 
states of the system. The ratio < /^(l, 2)// M (l, 1) < 1 
measures the site correlation of quanta: it is small when 
quanta are site-correlated (when there are many quanta 
on one junction there are almost none on the other one) 
and close to one otherwise. 

In Figl5]we show the nearest neighbor energy spacing 
(splitting) and the correlation function of the eigenstates. 
For this, and all the rest, we used I c = 13.3 /iA, Cj = 4.3 
pF, and ( — 0.1, which are typical values in experiments. 
We see that in the central part of the spectrum the en- 
ergy splitting becomes small in comparison to the aver- 
age. The corresponding pairs of eigenstates, which are 
tunneling pairs, are site correlated, and thus QBs. In 
these states many quanta are localized on one junction 
and the tunneling time of such an excitation from one 
junction to the other (given by the inverse energy split- 
ting between the eigenstates of the pair) can be expo- 
nentially large and depends sensitively on the number of 
quanta excited. 
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FIG. 2: (a) Energy splitting and (b) correlation function vs. 
energy for the two-junctions system (open circles, symmet- 
ric eigenstates; filled circles, antisymmetric eigenstates) . The 
labeled arrows mark the energy corresponding to the peak 
of the spectral intensity in Fig(3}b, d, and f (see text). The 
parameters are 7 = 0.945 and £ = 0.1 (22 levels per junction). 



The fact that the most site correlated eigenstates occur 
in the central part of the energy spectrum may be easily 
explained as follows: Let N be the highest excited state in 
a single junction, with a corresponding maximum energy 
AU (Fig|T|) . For two junctions the energy of the system 
with both junctions in the N-th state is 2AU, which 
roughly is the width of the full spectrum. Thus states 
of the form \N,0) and |0,iV) that have energy AU are 
located roughly in the middle. 

With the eigenvalues and eigenstates we compute the 
time evolution of different initially localized excitations, 
and the expectation value of the number of quanta at 
each junction (rii)(t) = (\&(i)|nj|\&(£)). Results are shown 
in FigUJi, c, and e. Also we compute the spectral in- 
tensity 1° = Kx^l^o)! 2 ! which measures how strong the 
initial state |^o) overlaps with the eigenstates. Results 
are shown in Fig[3]-b, d, and f, where we can see a peak 
in each case, which corresponds to the arrows in FigJ^fb. 
We can see that the initial state |^o) = 1 0, 5) overlaps 
with site correlated eigenstates with an energy splitting 
between them being relatively large and hence the tun- 
neling time of the initially localized excitation is short. 
For the case |\I>o) = |0, 19) QBs are excited: The excita- 
tion overlaps strongly with tunneling pairs of eigenstates 
in the central part of the spectrum, which are site corre- 
lated and nearly degenerate. The tunneling time of such 
an excitation is very long, and thus keeps the quanta lo- 
calized on their initial excitation site for corresponding 
times. Finally the initial state |^o) = |9, 19) overlaps 
with weakly site correlated eigenstates with large energy 
splitting. Hence the tunneling time is short. 

We tested whether a (coherent or incoherent) spread- 
ing of the initial state over a suitable energy window af- 
fects the results discussed above. For instance instead of 
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FIG. 3: Time evolution of expectation values of the num- 
ber of quanta at each junction for different initial excita- 
tions and their corresponding spectral intensity, (a) and 
(b): |*o) = |0,5); (c) and (d): |* ) = |0, 19); (e) and (f): 
\&o) = |9, 19). Open circles, symmetric eigenstates; filled cir- 
cles, antisymmetric eigenstates. The energies of the peaks in 
the spectral intensity are marked by labeled arrows in Fig[2]-b 
(see text). The parameters are 7 = 0.945 and £ = 0.1 (22 
levels per junction) . 



using the basis state |0, 19) as the initial state (Fig[3]-c,d), 
we superposed the basis states |0,20), |0, 19), |0, 18), and 
|0, 17). We found that the results qualitatively do not 
change. 

The experimental observation of QBs may be possible 
using the scheme of McDermott et al for simultaneous 
state measurement of coupled Josephson phase qubitsj^], 
where by applying current pulses the time evolution of 
the occupation probabilities in the qubits is measured. 
By applying a first microwave pulse on one of the junc- 
tions we excite it into a high energy single-junction state 
with energy ei and leave the other one in the ground 
state. In this way we have an initial state similar to 
the ones shown in Fig|3J After a variable time period 
we apply simultaneous current pulses to the junctions to 
lower their energy barriers AU and enhance the probabil- 
ity of tunneling outside the potential well. Then we test 
which junction switches to the resistive state (detected 
by a measurable voltage accross it). By repeating the 
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measuring many times we obtain the populations in the 
junctions as a function of the time period between the 
initial pulse and the simultaneous measuring pulses. 

Let us discuss the so far neglected quantum escape. 
For that we computed T escape by using the semiclassical 
formula [33| 

T eslape(£) = exp I - - J p((f)dip \ , (9) 

where a and b are the turning points of the classical 
motion in the reversed potential at U((p) — e, p(<p) — 
\/2\U (ip) — i], and Lo{e)/2-K is the frequency of the os- 
cillations inside the initial well. In table |l] we show 
the escape time from different metastable states, and 
we compare it with the tunneling time Ttunnei of an ini- 
tial excitation |W(0)) = |0, 1) between the two junctions, 
estimated from the energy splitting of the (symmetric- 
antisymmetric) pair of eigenstates with the largest over- 
lap with the initial excitation. We see that for Z = 19, 
where we excite QBs, the escape time is long enough 
for observing at least one tunneling exchange between 
the two junctions before escaping to the resistive state. 
Note that the cases I = 18 and 17 also excite QBs which 
would show even more tunneling exchanges before escap- 
ing. The case I — 16 does not excite QBs but eigenstates 
that, though having small energy splitting, do not show 
strong site correlation of quanta as in the previous cases. 
From these results we expect that escaping to the resis- 
tive state will not prevent from the experimental obser- 
vation of QB excitations. 

Another phenomenon that was not taken into account 
in our quantum model is decoherence. To be able to ob- 
serve tunneling between the junctions the coherence time 
has to be longer than the shortest tunneling time between 
the junctions, which is on the order of 1 ns in the cases 
shown in Fig|5]-a and e. In the experiment shown in [24| 
using a few levels per junction they obtained a coherence 
time on the same order. However, in the experiment in 
@ the coherence time was about 25 ns, and more recently 
in [tJ the coherence time was approximately 80 ns. We 
expect that further improvements in experiments Q will 
give us longer coherence times. 

Note that the above coherence times are shorter than 
the tunneling times of QB excitations (see table |T|, hence 
decoherence is an effect that can not be ignored if one 
wants to do a more realistic quantum description of the 
system. When exciting a JJ to a high-energy state, re- 
laxation (over dephasing) is usually the main source of 
decoherence. We can make a crude estimation of the cor- 
responding relaxation time T\ by using T\ ~ hQ/ei (Q 
is the quality factor of the junctions), which holds for a 
harmonic potential [34], [35|. For I = 19, 18 and 17, ei[h 
is around 150 GHz (see FigO-b). For the JJs used in y|, 
Q is between 500 and 1000, which leads to a relaxation 
time between 3 ns and 6 ns. It is much smaller than 
the tunneling time of the QB excitations, therefore one 
would expect to see instead of tunneling, a freezing of the 



TABLE I: Escape times for metastable states in a single JJ 
Tescape estimated by formula @, and tunneling time of the ini- 
tial excitation|*I'(0)) = |0, 1) between the two junctions r ttlnne ; 
estimated from energy splittings. 
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QBs on one of the junctions before they decohere due to 
relaxation. 

One could obtain more feasible results by increasing 
the bias current in such a way that there are less energy 
levels in the junctions. With this, exciting a QB would 
need less energy, and the relaxation time becomes longer. 
The tunneling time of that QB excitation is shorter, and 
might be even shorter than the relaxation time, allowing 
one to observe tunneling before relaxation. This possi- 
bility, and the inclusion of decoherence in our model, are 
issues that will be addressed in a future work. 

In summary, we have studied the classical and quan- 
tum dynamics of high-energy localized excitations in a 
system of two capacitively coupled JJs. In the classical 
case the equations of motion admit time periodic local- 
ized excitations (discrete breathers) which can be numer- 
ically computed. For the quantum case we showed that 
excitation of one of the junctions to a high level leaving 
the another junction in the ground state lead to a QB 
with long tunneling time. This is possible because the ex- 
citation overlaps strongly with tunneling-pair eigenstates 
which live in the central part of the energy spectrum and 
localize energy on one of the junctions. This result would 
not qualitatively change if we excite a (coherent or inco- 
herent) superposition of several product basis states in- 
stead of only one. We showed that with the available 
techniques for manipulating JJ qubits the experimental 
observation of QB excitations is possible. Escaping to 
the resistive state of the junctions (which together with 
decoherence was not taken into account in our quantum 
model) would not prevent us from doing that, and we 
expect that improvements in preparation (higher quality 
factors) and isolation techniques of JJ will lead to long 
enough coherence times, such that the phenomena we de- 
scribed in this work will be clearly observed. That would 
ultimately pave the way of a controlled stirring of quanta 
on networks of JJs. 
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